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Abstract. We analyze Ha observations of the recurrent nova T CrB obtained during the last decade. For the first time the 
Hot emission profile is analyzed after subtraction of the red giant contribution. Based on our new radial velocity measurements 
of the Ho- emission line we estimate the component masses of T CrB. It is found that the hot component is most likely a 
massive white dwarf. We estimate the inclination and the component masses to be i =i 67°, M WD = 1.37 ± 0.13M o and 
Mrg — 1.12 ± 0.23M o , respectively. The radial velocity of the central dip in the Ha profile changes nearly in phase with that 
of the red giant's absorption lines. This suggests that the dip is most likely produced by absorption in the giant's wind. 
Our observations cover an interval when the Ha and the (/-band flux vary by a factor of ~ 6, while the variability in B and V 
is much smaller. Based on our observations, and archival ultraviolet and optical data we show that the optical, ultraviolet and 
Ho- fluxes strongly correlate. We argue that the presence of an accretion disc can account for most of the observed properties of 
TCrB. 
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1. Introduction 

T Coronae Borealis is a well known recurrent nova with an 
orbital period of 227 d 57 (Kenyon & Garcia 1986, Fekel et al. 
I2000i . The system contains a M4.5 III giant (Murset & Schmid 
1999 1, which fills its Roche lobe and transfers mass toward a 
hot component (HC). The optical spectrum of T CrB is domi- 
nated by the red giant (RG) with superimposed Balmer emis- 
sion lines, which are thought to originate in an accretion disc 
(AD) around the HC (Selvelli et al. 1992; hereafter lSCGl . The 
high excitation Heii/14686 emission line is sometimes also 
visible (Iijima 119901 . The Ho- emission line is well visible 
all the time and shows variability with the orbital period as 
well as large variations on a time scale of years (Anupama & 
Prabhu 1991 Zamanov & Marti 2001 1. The photometric vari- 
ability in the long wavelength range is dominated by the ellip- 
soidal variability of the RG (Shahbaz et al. 1997 Belczynski 
& Mikolajewska 1998 , Yudin & Munari 1993 1. At short wave- 
lengths T CrB shows strong flickering activity, with the same 
parameters as those of normal cataclysmic variables, in spite 
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of the large difference in the size of the AD (Zamanov & 
Bruch 1998 1. Apart from the recurrent nova outbursts, T CrB 
shows other periods of enhanced luminosity. The amplitude 
of these variations increases toward the shorter wavelengths 
and reaches ~ 2 mag in the f/-band (Zamanov & Zamanova 
119971 Hric et al. I1998i . The equivalent width of the line Ha 
(EW(Ha)) also greatly increases and reaches ~ 30 - 35 A 
(Anupama & Prabhu IT99T1 Zamanov & Marti |200TT |. This ac- 
tivity is not well understood. 

For more than 30 years the HC in T CrB has been consid- 
ered as a main sequence star (Kraft 19581 Kenyon & Garcia 
1986 1. Investigations over the last 10 years, however, suggest 
that the accreting object must be a white dwarf (WD). First, 
ISCGI analvzed long-term ultraviolet (UV) spectral observations 
of T CrB and concluded that the HC is a WD. Belczynski & 
Mikolajewska ( 1998 1 analyzed the ellipsoidal variability of the 
RG and suggested the presence of a WD. In a study of recur- 
rent novae with giant secondaries, Anupama & Mikolajewska 
( 1999) also found that the HCs in these systems are WDs em- 
bedded in the dense wind of the RG. Finally, Hachisu & Kato 
( 2001 ) successfully modeled recurrent novae outbursts as ther- 
monuclear runaway on a massive WD. However, no system 
parameters determination based on the radial velocities (RVs) 
of the two components has been published since Kraft ( 1958 1. 
Recently, Hric et al. (119981 re-analyzed the RVs of Kraft J1958i 
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Fig. 1. Ha emission line profiles before (left panels) and after (right panels) the subtraction of the giant spectrum. The first number 
in each of the right panels shows the date of observations in YYMMDD format, and the second the orbital phase. Note that these 



phases were calculated with the orbital period of of Fekel et al. (2000 1, but the zero phase refers to the inferior conjunction of the 
RG. 



and concluded that the hot component in T CrB is a WD. 
However their analysis is still based on Kraft's seven measure- 
ments of the HC RV. 

In this paper we present an analysis of spectral Ha and pho- 
toelectric UBV observations of T CrB obtained during the last 
decade. This is the first investigation of this object where the 
RG contribution is subtracted and the "cleaned" Ha profile is 
measured. Based on our Ha radial velocity measurements we 
obtain a new orbital solution. We also discuss the long-term 
variability of T CrB. 

2. Observations and spectra processing 

The star T CrB was observed during 64 nights between 
February 1993 and September 2001 with the Coude spec- 



trograph of the 2.0 m RCC telescope of the Bulgarian 
National Astronomical Observatory "Rozhen". The spectra 
cover ~ 100 A or ~ 200 A around Ha, with resolution of 
~ 0.2 A pixel - . In some nights we have also obtained spectra 
covering the Heii/14686 emission line. Usually 2 or 3 expo- 
sures, of 15-20 minutes each, were taken to improve the sensi- 
tivity and to avoid problems with cosmic ray hits. The spectra 
were reduced in the standard way including bias removal, flat- 
field correction, wavelength calibration and correction for the 
Earth's motion. The spectra obtained within each observational 
season (with a typical duration of 2-4 days) were averaged to 
improve the signal-to-noise ratio. 38 mean Ha spectra were ob- 
tained. They are plotted in the left panels of Fig.^ We have 
also obtained new photoelectric UBV observations with the 60- 
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Table 1. The equivalent width EW and flux of the Ha line, 
its equivalent width after the giant spectrum subtraction EW(- 
gM), the factors /, and the radial velocity data of the HC, RG 
and central dip. The accuracy of the measurement is ~ 10% 
in EW (mainly due to the uncertainty of the continuum place- 
ment), ~ 20 - 40% in the flux (continuum placement and un- 
certainty of the continuum flux around Ha) and EW(-gM) (un- 
certainty of /) and ~ 10 - 15 km s _1 in the radial velocities. 
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cm telescope at Rozhen Observatory (Zamanov et al. 2004, in 
preparation). 

The main source of the continuum around Ha is the RG 
and the contribution of the HC is rather small, of the order of 
5-15%. Because of this and since the Ha emission line is usu- 
ally weak, the profile is strongly affected by the absorption lines 
of the RG. Therefore, to measure the RVs of the Ha emission 
line accurately one should first subtract the spectrum of the gi- 
ant. To do this we used the following procedure. The M4III 
star HD 134807 (Buscombe 1999 1 was observed with the same 
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Orbital phase 

Fig. 2. The RVs of the giant, HC and the dip folded with the 



ephemeris of Fekel et al. (2000 1. The best fits are shown with 
dotted lines. 



instrumental setup as T CrB and served as a template. The tem- 
plate was continuum-normalized and shifted to zero RV. Once 
Ha was masked, the RVs of the giant were obtained by con- 
volving the T CrB spectra with the template (Tonry & Davis 
IT9791 . The template was then shifted to the RV of each spec- 
trum, multiplied by a number / between and 1 and subtracted 
from the T CrB spectra. / was chosen to minimize the scat- 
ter of the residuals. The width of the absorption lines in the 
spectra of both T CrB and the template is dominated by the 
spectral resolution rather than the rotational broadening. As a 
result the absorption lines have almost equal widths; however, 
in some cases the line widths are different. In those cases the 
average absorption line widths in the template and T CrB spec- 
tra were estimated by an autocorrelation function. Before the 
subtraction, the template or the spectrum was broadened by 
convolving it with a Gaussian to equalize the line widths. The 
"cleaned" Ha emission line profiles and the factors / used are 
shown in the right panels of Fig.^and in Tabled respectively. 

3. Radial velocities 

The radial velocities of the Ha emission line were measured 
by fitting the line wings with a Gaussian function. We fitted 
only the outer parts of the line which are not contaminated 
by the central asymmetry/dip. Figure^ shows that even af- 
ter subtraction of the giant's spectrum the Ha line is double- 
peaked. When possible, the RV of the dip was also mea- 
sured. A simple fitting with two Gaussian functions was used 
- one broad with a positive amplitude for the emission line 
and another narrow with a negative amplitude for the dip. 
Figure|2] shows the RVs of the giant, Ha emission line and 
the dip folded with the ephemeris of Fekel et al. ( 2000 1 To — 
/D2447918.62 + 227 d 5687£. The RVs were fitted with the 
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Table 2. The fit parameters to the Ha RV data. 





red giant 


hot component 


central dip 


7 [km s 1 ] 


-28.7 (1.4) 


-36.7 (1.6) 


-36.1 (2.3) 


K [kms- 1 ] 


23.7 (1.8) 


19.5 (2.1) 


15.2 (3.0) 


00 


-0.009 (0.013) 


0.563 (0.019) 


-0.10(0.03) 


cr [kms -1 ] 


8.2 


9.0 


10.0 



function V, = y+/Tcos[2^(0-0o)] and the obtained parameters 
are given in Table|2] 

The RV of the giant follows the ephemeris of Fekel et al. 
(2000) very well. The RV of the Ha wings changes in anti- 
phase with that of the giant. This suggests that the location of 
the emission is on the line connecting the two stars, on the 
side of the mass center which is opposite to the giant. Most 
likely, the region where Ho- emission line is formed is centered 
close to the HC. The observed ellipsoidal variability (Shahbaz 
et al l 19971 Yudin & Munari 1993 Belczynski & Mikolajewska 
1998 1 suggests that the giant fills its Roche lobe, which pro- 
vides a possibility for mass flow through the inner Lagrangian 
point and an AD formation. In this case the most plausible 
assumption is that the bulk of the Ha flux is emitted by the 
outer part of the AD. There is also other indirect evidence 
for the presence of an AD in this system ( SCG Anupama & 
Mikolaiewska fT999l Hachusu & Kato l200H . Since the AD or- 
bits around the mass center of the system together with the HC, 
the RVs of the emission line wings are likely to trace the mo- 
tion of the HC. We will then use the radial velocity curves to 
estimate the component masses. Note, however, the small dif- 
ference in the y velocity of the HC and RG, and that <po of the 
HC RVs is not exactly 0.5. This is likely to be a result of the 
presence of additional emission components which most prob- 
ably arise in the ionized parts of the giant's wind. Such emis- 
sion components can be easily seen in Fig.[0 (the panels for 
960608, 970618 and 980811). 

If the double-peaked shape of Ha is a result of emission 
from an AD seen under moderate-to-high inclination, then the 
RV of the central dip should be in phase with that of the line 
wings. Figure|2]shows that the RV of the dip is modulated with 
the orbital period but follows more closely that of the giant, 
which argues against the AD origin. The central dip in the 
emission profiles of the symbiotic stars is thought to originate 
mainly in the giant's wind. In the symbiotic star EG And the RV 
of the dip also changes in phase with that of the giant (Munari 
119931 . Thus, the central dip in the Ha emission profile of T CrB 
originates most probably in the giant's wind. 

4. Orbital solution 

Having the RV amplitudes of the two components we can es- 
timate their masses. The velocity amplitude of the giant in 
T CrB was accurately determined by Kenyon & Garcia ( 1986 1 
- ^rg = 23.32 + 0.16 kms"' and was also recently revised 
to K RC = 23.89 + 0.17 kms" 1 by Fekel et al. J2U001 ). Our 
giant's velocity data are in full agreement with the results of 
these authors. But since our measurements are noisier, they 
lead to larger uncertainty of K R q- Thus, for determination of 
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Fig. 3. The masses of the components of T CrB as a function of 
the inclination. The dotted curve shows the most likely position 
of the HC and the i - mass space. 

the system's parameters we prefer to use K R q of Fekel et al. 
(2000). Assuming a zero eccentricity we obtain the following 
parameters: Mhc sin 3 ; = 1.06 + 0.10M o and Mrg sin 3 / 
0.87 ±0.18M o . 

The masses of the two components as a function of the 
system inclination are shown in Fig.[3] Our HC velocities al- 
low solutions with Mhc below the Chandrasekhar limit when 
the inclination is i > 65°. To obtain a more precise estima- 
tion of the masses a tighter constraint on the inclination is 
needed. The HC contributes little in the visible wavelengths 
but it gives almost the whole UV flux. Therefore, if the HC 
was eclipsed this would be best observed in the UV. The UV 
data, however, does not show any hint of eclipse ( SCG I and it 
can provide an upper limit. We assume that the giant fills its 
Roche lobe. We have then numerically calculated the projec- 
tion of the RG in the orbital plane at zero phase for several 
values of i between 60° and 80°, and component mass ratio 
q = M RG /M hc = ^hc/^rg = 0.82 + 0.10. We found that 
the center of the HC is not eclipsed when i < 70°. The upper 
limit depends also on the size of the UV emitting region. Its 
size, however, should be very small compared to the radius of 
the HC's Roche lobe because the most likely UV emitting re- 
gion is the WD and the inner part of the AD (SCG I. Thus, we 
conclude that the upper limit of the inclination is ~ 70°. The 
dotted lines in Fig.|3]mark the limits of the inclination and +lcr 
of the HC mass. Recently, Hachisu & Kato (2001 1 have mod- 
eled the outburst light curve of T CrB and obtained the WD 
mass Mwd - 1-37 + 0.01M o . In terms of our orbital solution 
this corresponds to i ^ 67° and for this inclination we have 
M WD a 1.37 + 0.1 3M and Mrg ~ U2 + O.23M . 

5. Long-term variability 

In Fig.|4] we show our observations of T CrB along with 
those published. Figure^ shows our EW(Ha) measurements 
(dots) and those of Anupama & Prabhu A199U (open cir- 
cles). Figures|4j> and c show the total UV flux between 1250 
and 3200 A (F uv ; taken from lSCGl and EW(Ha) measured 
from our spectra after the subtraction of the RG contribution 
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Fig. 4. a) Ha equivalent width; b) total UV flux between 1250 and 3200 A; c) Hot equivalent width after the subtraction of the 
RG spectrum; d) U magnitudes, e) B magnitudes and f) V magnitudes. See text for the sources of the different quantities. 



(hereafter EW(Ha-gM)), respectively. In the photometry pan- 
els (Figs.@J/ - f) the dots show the standard UBV magni- 
tudes (Raikova & Antov 119861 Lines et al. 119881 Hric et al. 
1998 Sherrington & Jameson 1983 Zamanov & Zamanova 
119971 our new data) and the open circles the instrumental ones 
(Oskanian [T9^3l Luthardt [1 992 1. The data of Luthardt were 
presented only in graphical form as a difference between T CrB 
and the comparison star which was not specified. We estimated 
the magnitudes of T CrB from the figure of Luthardt ( 19Sj2) in 
the following way. On the computer screen we pointed with the 
cursor on each data point and read its coordinates. These "com- 
puter screen" coordinates were converted to Julian Dates and 
magnitudes by applying the necessary transformations. Then, 
the data were shifted to the scale of the standard magnitudes by 
using the intervals when they overlap with the U BV observa- 
tions. We have also obtained an estimation of the Ha flux by 
using the measured EW(Ha) and the published /?-band mag- 
nitudes. We assume that the R-band flux represents the contin- 
uum flux at the position of Ha and that the variations in R are 
entirely due to the ellipsoidal variability. We fitted the R data of 
Mikolajewski et al. (119971 and Hric et al. d!998i with a func- 
tion consisting of three sines with periods of P or b, P rt>/2 and 
^orb/3. With the coefficients obtained we have calculated the R 
magnitudes at the moments of the spectral observations. Then, 
Bessell's ( 1979 1 calibration of the /?-band was used to calculate 



the Ha flux. Because the Ha flux variations repeat very closely 
those of EW(Ha), the Ha flux is given only in Tabled (only 
for our spectra). We note that the assumption that the /?-band 
flux represents the flux in the continuum around Ha introduces 
the largest uncertainty in the calculated Ha flux. The reason is 
that the R filter covers several strong molecular bands, and is 
placed where the flux in the RG spectrum increases strongly 
toward the long wavelengths. 

As seen from Fig.0]our observations cover an interval when 
both the photometric data and EW(Ha) vary greatly. The long- 
term variability has amplitude ~ 2.0, 1.0 and 0.3 mag in U, B 
and V bands, respectively. EW(Ha) increased from ~ 5 A to 
~ 30 A. EW(Ha-gM) also increased, but its maximum is de- 
layed with respect to the maximum of EW(Ha) and appears at 
JD2450950. It is seen from Figs.@]af and e that the maximum 
of EW(Ha-gM) coincides with the rapid decrease of the U and 
Z?-band fluxes. This might explain the observed peculiarity. It 
is also seen from Fig.|4]that most of the UV observations have 
been obtained when T CrB was in a high state. It is quite evi- 
dent, however, that during the observations around JD2444000 
and JD2447700 T CrB was in a low state. The average increase 
of Fuv is ~ 3.3 and the corresponding increase of the [/-band 
flux is ~ 2.6 (~ 1 mag). 

Anupama & Mikolajewska ( 1999 1 noticed that the increase 
of EW(Ha) in the late 1980s correlated with the increase of UV 
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and U fluxes. Our new observations also show similar behav- 
ior. To verify the strong visual impression for correlation of the 
data shown in Fig.0]we performed a correlation analysis. The 
small number of UV observations required interpolation when 
calculating the correlations involving F\jy, For the correlation 
between Ha and U -band fluxes we used yearly averaged val- 
ues. (EW(Ha-gM) was not used here). The computed Pearson's 
correlation coefficient and Spearman's rank correlation coeffi- 
cient are in all cases larger than 0.7 with the significance of the 
correlation always higher than 99%. The correlation analysis 
confirmed that all the data in Fig.|4]are indeed strongly corre- 
lated. As an example we show Ha vs. t/-band flux in Fig.|5] 
Our data and those of Anupama & Prabhu ( 1991 ) follow differ- 
ent tracks, but the correlation is apparent. 

The correlated variations in the continuum can be also 
seen in Fig. [6] where we present the UV-optical-IR spectral 
energy distribution (SED) of T CrB in low and high states. 
The points are a compilation (average) from the photometric 
data of Shahbaz et al. fl997l . Kamath & Ashok ( fl^99l . Yudin 
& Munari JT993b . Munari et al. (fT992t . Hric et al. JT9981 . 
Mikolajewski et al. J1997i . Zamanov & Zamanova ( 1997) and 
our measurements. A representative high-state UV spectrum of 
T CrB was obtained by averaging the low-resolution spectra 
taken by the International Ultraviolet Explorer {WE) between 
JD2444300 and JD2447300 when the system was in high state. 
The rest of the IUE spectra were used to obtain the average 
low-state UV spectrum. All fluxes were corrected for interstel- 
lar reddening of Eb-v = 0.15 (Cassatella et al. 1982) using 
O'Donnell's ( 1994) interstellar extinction law and Ry = 3.1. 

The data in Figs.|4j> and[6]show that the bolometric lumi- 
nosity of the hot component varies greatly on time scale of 
about 200 days, which is not a typical feature of objects whose 
activity is powered by thermonuclear runaways. Moreover, 
ISCGI demonstrated that the bolometric luminosity of the hot 
component is low (~ 100L o ) compared to some other symbi- 
otic systems (e.g., AG Peg, AG Dra, Z And ) whose activity 
is supposed to be due to hydrogen burning and that even in 
the high state most of the hot component's light is emitted in 
the UV. All these characteristics suggest that the system is ac- 
cretion powered and the hot component is most probably the 



innermost part of an optically thick accretion disc and/or the 
boundary layer between the disc and the WD's surface as was 
already proposed by SCG This idea supports our previous as- 
sumption based on the behavior of the Ha- radial velocity (Sect. 
3) that the bulk of the energy of this line is probably emitted by 
the outer part of an AD. 

According to the AD theory (e.g., Warner 1995), there ex- 
ists a critical value M cr of the mass accretion rate M, so that 
if M > M cr the disc is in a steady state, being hot and opti- 
cally thick with effective temperature r e ffj,ot x M'' 4 M^r 3 ' 4 
(TW.hot > 10 4 K). Here r is the distance to the WD. If M < M cr 
the disc is cool and optically thin with r e ff, C ooi — 5 - 6 x 10 3 K. 
With the parameters of T CrB M ~ lO^Moyr" 1 l ISCGT ) and 
Mwd - 1-37, only the inner part of the disc with radius ~ IR Q 
is in the hot state, while the rest of the disc is cool. In the case 
of accretion at a high rate onto a very massive WD, as in T CrB, 
the temperature of the innermost part of the disc should be 
~ 1 - 2 x 10 5 K. It will thus emit mostly in the UV and can 
be the HC. In this case variations of M can easily account for 
its observed variability. 

The symbiotic systems containing a hot compact object 
change their energy distribution appreciably during active 
phase because of the variation of the temperature of the com- 
pact object. Figure|6] shows that the UV energy distribution of 
T CrB tends not to vary during its transition from low to high 
states. ISCGI notice a correlation of the slope of the contin- 
uum with the total UV flux. However, using the values given 
in their Table 1, we found that the correlation, if any, is very 
weak. Thus, the modeling of the continuum requires nearly 
the same temperature of the hot component in low and high 
states. We tried to model the continuum of T CrB assuming 
the presence of an optically thick AD. The IR region was fit- 
ted with synthetic fluxes of a model spectrum of a red giant 
with r eff = 3500 K, M = 1M (P. Hauschildt, private commu- 
nication) and R = 15R Q , filling its Roche lobe. The UV and 
visual regions were fitted with an emission of an AD and gas 
environment. We adopted an electron temperature of 10000K 
according to lSCGI Thus, the only free parameters are the mass 
accretion rate and the relative scaling between the three com- 
ponents. 

To determine the nebular continuum, the ionization state 
of helium has to be known. The dominant state of ionization 
can be estimated from the ratio of the emission measures of 
the neutral and ionized helium based on the fluxes of some 
lines of Hei and the Heii/14686 line. Since the electron tem- 
perature of 10000 K is comparatively low we assumed that 
He i lines are determined only by recombination. Our obser- 
vations of the He n /14686 line show that most of the time this 
line was absent and only occasionally appeared very weak. The 
line was also not presented in the spectrum in July 1986, as 
can be seen in Fig. 2 of ISCGI According to Iijima (1990) in 
the period 1981-1987 T CrB has had prominent emission lines 
of He i whereas its He n /14686 line was only temporarily ob- 
served. In the spectra taken in Apr-May 1997 during the re- 
cent high state, Anupama & Mikolajewska ( 119991 ) detected the 
He ii /14686 line, but it was by a factor of 2 weaker than the 
He i lines. At an electron temperature of 10000 K and elec- 
tron density of 10 6 cm 3 if the number of He + and He ++ ions 
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Fig. 6. UV-optical-IR energy distribution of T CrB. The representative UV spectra are related to its low and high states. The 
triangles and the filled dots indicate the photometric fluxes at low and high states, respectively. Lines a and b show the model 
spectrum of the system in high and low state, respectively. Line c connects the synthetic fluxes of a red giant. Lines d and e 
indicate the continia of the accretion disc and the nebula in high state. Lines f and g show the same components for the low state. 



is the same, the intensity of the Heii/14686 line exceeds the 
intensities of the Hei lines by one order of magnitude (e.g. 
Pottasch 1984i. This was probably not the case during the pe- 
riod of Iijima's observations. Thus, we accepted the dominant 
helium ionization state He n in all occasions and assumed that 
the nebular continuum is mostly continuum emission of hydro- 
gen and neutral helium. We adopted a helium abundance of 0. 1 
(Nussbaumer & Vogel 1989) and distance to the system of 1.3 
kpc (Patterson[T$84). 

The result of the modeling is shown in Fig. [6] To deter- 
mine the level of the UV continuum of T CrB is a complicated 
task - the spectrum is dominated by humps caused by blended 
emission lines. We fitted the bottoms of the troughs between 
those humps. The largest discrepancies between the fit and the 
spectrum are in the short- wave length region and around 2200 A 
where the sensitivity of WE is low. It is seen, however, that the 
shape of the fit to the UV spectrum in the two states is nearly 
the same. The observed optical fluxes in both states are sys- 
tematically higher than the fit. This indicates the presence of 
component(s) which are not taken onto account in our model. 
This additional light might be due to density and temperature 
inhomogeneities in the outer disc or/and in the other nebular 
formations in the system. 

The modeling gave an accretion rate of 4 x 10~ 9 M yr -1 for 
the low state and 3 x 10~ 8 M Q yr~' for the high state in agree- 
ment with lSCGl Note also that at these accretion rates the disc 
temperature exceeds 10 5 K. At such temperatures the scope of 
the UV continuum is nearly constant. 

We obtained appreciable emission measures of 2.7 x 
10 59 cirT 3 and 7.4x 10 59 crrr 3 for the low and high states, deter- 
mined from strong nebular UV and optical continuum. On the 
other hand the Balmer emission lines are weak (Ha exceeds the 
level of the continuum by a factor of not greater than 5-6). The 



symbiotic systems with an extended surrounding nebula (AG 
Dra, Z And) have prominent Balmer emission lines exceeding 
the continuum by a factor of up to 40-50, together with their 
intensive nebular continuous spectrum. In the case of T CrB 
the weakness of the Balmer emission lines can be understood 
if they are assumed to be formed in a gas medium with a great 
optical depth, which absorbs the dominant part of the photons. 
The optical depth is great when the column density of the emit- 
ting gas is high. Anupama & Mikolajewska ( 1999 1 came to the 
conclusion that the hot component of the symbiotic recurrent 
novae is most probably a white dwarf + accretion disc embed- 
ded in an optically thick envelope formed by the wind of the M 
giant secondary. This model provides high column density and 
contains an AD as well. So, the weak Balmer lines of T CrB 
can be qualitatively explained by the assumption of an origin 
in a binary system containing an AD. The investigation of the 
RV of the Ha line showed that the bulk of its energy is emitted 
by an area around the hot component - probably the outer part 
of an AD. This model is supported also by our suggestion for 
the explanation of the weak Balmer emission lines of T CrB. 
The modeling of T CrB's SED suggests that the UV spectrum 
appears flat mostly because of the strong nebular contribution. 
However, it is also possible that this is partly a result of seeing 
the AD at a relatively high inclination as suggested bv lSCGI 

The evidence for an AD and high M provides good con- 
ditions for formation of intensive high-speed wind emerg- 
ing from the innermost hot parts of the disc (see Mauche & 
Raymond 1997 for a discussion of disc wind in the case of cat- 
aclysmic variables). Consequently the nebular spectrum of the 
system may appear also in the winds of the disc and the cool 
giant in addition of the disc itself. Shore & Aufdenberg ( 1993 1 
showed that the giant's wind is variable. Further evidence for 
this can be found in the variations of the central dip of the Ha 
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emission line - most of the time the dip is well pronounced, but 
sometimes it is highly variable and is almost absent (Fig.0. 

If an AD is present in T CrB system, its activity may be a 
result of variation of the accretion rate M. There are two ways 
to increase M: (1) to increase the mass transfer rate M a from 
the giant and (2) to increase temporarily the mass flow through 
the AD. The former may be a result of the magnetic activity of 
the giant as suggested by Soker (2002). The later mechanism 
corresponds to the outbursts observed in the dwarf novae and is 
related to the so-called disc instability. In the context of symbi- 
otic stars the disc instability was studied by Duschl ( 1986 1, but 
assuming that the accreting object is an ~ 1M main-sequence 
star rather than a WD. Briefly, in disc the instability model 
M a is constant and less than M cr , i.e. the disc is in the cool 
state most of the time. Outburst is observed when M (locally in 
the disc) increases above M CI and the AD rapidly approaches 
the hot state (for more details see Warner 1995). In both cases 
the disc luminosity increases because of the dependence of the 
temperature on the mass accretion rate: r e ff,hot x M 1 ^ 4 - 

6. Summary and conclusions 

The radial velocities of the Ha emission line of the recur- 
rent nova T CrB were precisely measured after subtraction of 
the red giant spectrum, thus obtaining the unbiased Ha pro- 
file. Based on our 38 new radial velocity measurements we 
obtain a new solution for the system parameters. The solu- 
tion indicates that at a system inclination i > 65° the mass of 
the hot component is below the Chandrasekhar limit and thus 
the hot component in T CrB is most likely a massive white 
dwarf. The most likely values of the system inclination and the 
component masses are i 67°, Mwd - 1-37 + 0.1 3M and 
M RG ^ 1.12 + O.23M . 

The radial velocity of the central dip of the Ha emission 
line profile changes almost in phase with that of the red giant. 
This suggests that the central dip originates from the giant's 
wind. 

During our observations T CrB showed a period of activ- 
ity when the system brightness increased by ~ 2.0, 1.0 and 
0.3 mag in U, B and V bands, respectively. The Ha equivalent 
width also increased by a factor of ~ 6 reaching ~ 30 A and the 
line had a complex multi-component structure. The phasing of 
the modulation of Ha emission line radial velocity shows that 
the emitting region is gravitationally connected to the hot com- 
ponent of the system and most probably related to the outer part 
of an accretion disc. We argue that a number of other charac- 
teristics of T CrB are also in support of this supposition. The 
long-term variability is then assumed to be caused by varia- 
tions of the mass accretion rate, possibly due to changes of the 
mass transfer rate from the companion or disc instability. This 
is also supported by our success in modeling the UV-optical- 
IR SED of T CrB with a three-component model: an opti- 
cally thick accretion disc, a nebular continuum and a red giant. 
However, to construct a consistent picture and to understand 
the underlying physical processes in T CrB better, simultane- 
ous UV and optical spectral observations during future active 
phases are needed. Space observatories offering simultaneous 
optical/UV spectroscopic capabilities like Rronos (Peterson et 



al-HSD and WSO/UV (Wamsteker & Ponz 2001 1 would be a 
good solution. 
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